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The  recent  study  of  the  sensitivity  of  tropical-cyclone  intensification  to  the  surface 
drag  coefficient  in  a  three-dimensional  model  by  Montgomery  et  al.  is  extended  to 
include  a  wind-speed-dependent  drag  coefficient  and  one  of  four  boundary-layer 
parametrization  schemes:  the  bulk,  Blackadar,  MRF  and  Gayno-Seaman  schemes. 
The  schemes  are  slightly  modified  to  have  the  same  drag  coefficient  formulation 
and  the  same  constant  exchange  coefficients  for  sensible  heat  and  moisture.  Interest 
is  focussed  on  the  change  in  intensity  of  the  azimuthally-averaged  tangential  wind 
speed  and  change  in  the  low-level  vortex  structure  when  the  standard  value  of  the  drag 
coefficient  is  halved  or  doubled.  Changing  the  drag  coefficient  provides  insight  into 
unbalanced  effects  in  the  boundary  layer  and  their  impact  on  the  vortex  evolution 
and  structure.  The  changes  in  vortex  behaviour  with  changing  drag  coefficient  are 
qualitatively  similar  for  all  schemes,  the  maximum  intensification  occurring  for  a 
value  somewhere  near  the  standard  value  of  the  drag  coefficient.  The  interpretation 
given  to  explain  this  behaviour  underlines  the  intrinsically  unbalanced  nature  of  the 
boundary-layer  dynamics,  although,  for  reasons  discussed,  a  complete  theory  for 
the  behaviour  does  not  exist.  The  behaviour  found  is  at  odds  with  the  predictions  of 
Emanuel’s  (balance)  theory  for  the  maximum  intensity  of  a  tropical-cyclone,  which 
predicts  a  monotonic  decrease  in  intensity  with  the  drag  coefficient  if  the  enthalpy 
exchange  coefficient  is  held  fixed.  It  is  at  odds  also  with  a  recent  numerical  study  of 
the  maximum  intensity  by  Bryan  and  Rotunno. 

The  study  underscores  the  importance  of  boundary-layer  dynamics  in  models  for 
forecasting  tropical-cyclone  intensity  and  the  need  for  care  in  choosing  a  boundary- 
layer  scheme.  However,  it  is  not  yet  known  which  boundary-layer  formulation  is  the 
most  appropriate  for  this  purpose,  highlighting  the  need  for  a  concerted  research 
effort  in  this  direction. 

Key  Words:  hurricanes;  tropical  cyclones;  typhoons;  surface  drag  coefficient;  frictional  drag;  boundary  layer 

Received  16  June  2010;  Revised  9  August  2012;  Accepted  19  September  2012;  Published  online  in  Wiley  Online 
Library  22  January  2013 

Citation:  Smith  RK,  Montgomery  MT,  Thomsen  GL.  2014.  Sensitivity  of  tropical-cyclone  models  to 
the  surface  drag  coefficient  in  different  boundary-layer  schemes.  Q.  /.  R.  Meteorol.  Soc.  140:  792-804. 
D01:10.1002/qj.2057 


@  2013  Royal  Meteorological  Society 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

APR  2014  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2014  to  00-00-2014 

4.  TITLE  AND  SUBTITLE 

Sensitivity  of  tropical-cyclone  models  to  the  surface  drag  coefficient  in 
different  boundary-layer  schemes 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Department  of 

Meteorology, Monterey, CA, 93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

13 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Surface  Drag  Coefficient  and  Tropical  Cyclones 


793 


1.  Introduction 

A  tropical  cyclone  can  be  viewed  as  a  complex  heat  engine 
that  extracts  heat  energy  from  the  sea  surface,  mainly  in  the 
form  of  latent  heat,  and  radiates  sensible  heat  to  space  in 
the  upper  troposphere  (Kleinschmidt,  1951;  Emanuel,  2004; 
and  references  therein).  It  loses  mechanical  energy  primarily 
in  the  form  of  frictional  drag  at  the  ocean  surface  and 
does  work  in  creating  an  upper-level  anticyclone.  Angular 
momentum  associated  with  the  background  rotation  of  the 
earth  is  drawn  into  the  cyclone  by  an  overturning  circulation 
brought  about  by  the  collective  effect  of  buoyancy  in  deep 
rotating  clouds.  A  fraction  of  the  angular  momentum  drawn 
inwards  is  lost  because  of  the  frictional  torque  at  the  sea 
surface. 

Recent  research  using  three-dimensional  numerical 
simulations  has  shown  that  the  spin-up  of  intense  winds 
in  the  inner  core  of  the  cyclone  actually  occurs  within 
a  shallow  friction  layer,  or  boundary  layer,  near  the  sea 
surface,  which  is  approximately  500  m  deep  (Zhang  et  ah, 
2001;  Bui  et  al,  2009;  Smith  et  al.,  2009).  The  simulations 
indicate  that  the  spin-up  is  strongly  controlled  by  dynamical 
processes  in  the  boundary  layer  and  that  spin-up  in 
models  may  be  sensitive  not  only  to  the  representation 
of  surface  drag,  but  also  to  the  representation  of  vertical 
turbulent  mixing  of  momentum  in  the  boundary  layer.  To 
obtain  a  more  complete  understanding  of  these  sensitivities 
requires  idealized  modelling  studies  designed  to  isolate 
as  far  as  possible  the  basic  physical  processes  involved. 
In  this  spirit,  Montgomery  etal.  (2010)  re-examined  the 
sensitivity  of  tropical  cyclones  in  a  three-dimensional,  non¬ 
hydrostatic,  cloud-representing  numerical  model  to  the 
surface  exchange  coefficient  of  momentum  and  Smith  and 
Thomsen  (2010)  examined  the  sensitivity  to  the  choice  of 
the  boundary-layer  parametrization  with  the  same  exchange 
coefficients. 

Early  theoretical  and  numerical  studies  of  the  sensitivity 
to  the  surface  exchange  coefficients  in  axisymmetric 
models  found  that  the  intensity  decreases  markedly  with 
increasing  drag  coefficient  (Emanuel,  1995;  Craig  and  Gray, 
1996;  and  references  therein).  In  contrast,  Montgomery 
etal.  (2010)  found  that,  in  a  three-dimensional  model, 
the  intensification  rate  and  maximum  intensity  of  the 
vortex  increase  with  increasing  surface  drag  coefficient 
until  a  certain  threshold  value  is  attained,  and  then 
decrease.  Moreover,  whereas  a  vortex  intensifies  in  the 
Emanuel  and  Rotunno  (1987)  axisymmetric  model  when 
there  is  no  surface  drag,  no  system-scale  intensification 
occurs  in  the  three-dimensional  model  used,  despite 
persistent  sea-to-air  fluxes  of  moisture  that  maintain 
deep  convective  activity.  The  Montgomery  etal.  (2010) 
calculations  were  carried  out  using  the  Pennsylvania 
State  University-National  Center  for  Atmospheric  Research 
fifth-generation  Mesoscale  Model  (MM5)  with  the  bulk 
aerodynamic  boundary-layer  parametrization  and  constant 
(i.e.  wind-speed-independent)  surface  exchange  coefficients. 
Deep  convection  was  modelled  using  a  warm-rain  scheme 
and  radiation  was  parametrized  crudely  using  a  Newtonian 
relaxation  scheme  so  as  to  minimize  the  generation  of 
excessive  convective  instability  in  the  vortex  environment, 
where  the  grid  resolution  is  insufficient  to  represent 
convection.  The  new  findings  were  interpreted  using  recent 
insights  obtained  on  tropical-cyclone  intensification,  which 
highlight  the  intrinsically  unbalanced  dynamics  of  the 


tropical-cyclone  boundary  layer  and  the  three-dimensional 
nature  of  the  spin-up  process  (Nguyen  et  al.,  2008;  Bui  et  al, 
2009;  Smith  et  al,  2009). 

Smith  and  Thomsen  (2010)  carried  out  experiments  with 
a  similar  version  of  the  MM5  model  to  investigate  the  depen¬ 
dence  of  tropical-cyclone  intensification  and,  in  particular, 
the  structure  of  the  tropical-cyclone  boundary  layer,  on  the 
representation  of  the  boundary  layer  in  the  model.  Predic¬ 
tions  using  one  of  five  available  schemes  were  compared, 
not  only  between  themselves,  but  where  possible  with  recent 
observational  analyses  of  hurricane  boundary-layer  struc¬ 
ture.  These  schemes  included  the  bulk  aerodynamic  scheme, 
the  Blackadar  scheme,  the  Burk-Thompson  scheme,  the 
MRF  scheme  and  the  Gayno-Seaman  scheme.  The  same 
surface  exchange  coefficients,  modified  to  reflect  recent 
measurements,  were  used  in  aU  boundary-layer  schemes  to 
facilitate  a  proper  comparison  of  the  schemes.  The  study 
showed  that  there  is  a  significant  sensitivity  of  vortex  evolu¬ 
tion  to  the  particular  boundary-layer  scheme  used.  Sensitive 
aspects  include:  the  onset  time  of  rapid  intensification,  the 
low-level  wind  structure  in  the  inner-core  region,  and  the 
intensity  after  five  days  of  integration.  Some  of  these  differ¬ 
ences  were  traced  to  the  different  vertical  eddy  diffusivities 
determined  by  the  schemes.  In  particular,  the  MRF  scheme 
appears  to  be  overdiffusive  compared  with  the  others  in  the 
tropical  cyclone  context,  consistent  with  the  finding  of  Braun 
and  Tao  (2000)  and  also  with  the  few  available  observational 
estimates  at  high  wind  speeds  over  the  ocean  (Zhang  et  al, 
2011;  and  references  therein).  Nevertheless,  in  other  applica¬ 
tions  the  scheme  has  been  found  to  be  the  more  skilful  (e.g. 
Thomsen  and  Smith,  2008).  The  Smith  and  Thomsen  (2010) 
study  fell  short  of  being  able  to  advocate  the  use  of  a  par¬ 
ticular  scheme,  although  certain  shortcomings  of  individual 
schemes  were  identified  in  relation  to  their  ability  to  capture 
realistic  vertical  wind  profiles  and  surface  inflow  angles.  It 
was  pointed  out  that  the  current  inability  to  determine  ‘the 
optimum  scheme’  has  implications  for  the  predictability 
of  tropical-cyclone  intensification.  Nevertheless,  the  dif¬ 
ferences  between  schemes  provide  an  estimate  of  forecast 
uncertainty. 

Smith  and  Thomsen  (2010,  p  1682)  found  that  the 
schemes  with  the  smallest  vertical  eddy  diffusivity  have 
the  largest  inflow  speeds,  the  largest  tangential  winds,  and 
the  smallest  radius  at  which  they  occur.  Drawing  on  ideas 
emerging  from  calculations  with  a  slab  boundary  layer 
by  Smith  and  Vogl  (2008),  they  argued  that  the  effective 
frictional  stress  within  the  boundary  layer  increases  with 
decreasing  boundary-layer  depth,  because  then  the  surface 
stress  is  effectively  spread  over  a  shallower  layer.  Those 
calculations  showed  that  a  larger  effective  stress  leads  to 
a  larger  disruption  of  gradient-wind  balance  and  hence  to 
stronger  inflow  in  the  boundary  layer.  In  turn,  stronger 
inflow  brings  air  parcels  to  small  radii  more  rapidly, 
minimizing  the  loss  of  absolute  angular  momentum  on 
account  of  the  frictional  torque,  even  though  the  torque 
itself  is  increased.  In  this  way,  the  maximum  tangential 
wind  is  increased  relative  to  the  calculations  with  larger 
vertical  diffusivity.  The  finding  of  Smith  and  Thomsen 
concerning  the  effect  of  the  vertical  eddy  diffusivity  on  the 
maximum  azimuthally-averaged  tangential  winds  appears 
to  conflict  with  the  conclusion  of  Bryan  and  Rotunno 
(2009,  p  1777)  regarding  intensity  in  an  axisymmetric  model 
that  ‘there  is  a  strong  sensitivity  to  Zh>  but  essentially  no 
sensitivity  to  ly.  They  use  the  horizontal  and  vertical  mixing 
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lengths,  Zh  and  ly,  to  assess  the  dependence  of  cyclone 
intensity  on  the  subgrid-scale  mixing  formulation.  One 
could  argue  that  a  lack  of  sensitivity  to  ly  implies  a  lack 
of  sensitivity  to  the  vertical  eddy  diffusivity.  However,  one 
cannot  make  a  one-to-one  correspondence  between  the 
mixing  length  and  the  eddy  diffusivity  on  account  of  the 
assumed  dependence  of  the  diffusivity  not  only  on  the 
mixing  length,  but  also  on  the  local  flow  deformation  and 
static  stability. 

The  above  discussion  highlights  the  need  to  further 
understand  the  influence  of  the  boundary  layer  on  vortex 
intensification.  Given  the  efforts  that  have  been  devoted 
to  including  wind-wave  coupling  for  the  intensification 
problem  (Chen  et  al,  2007),  it  is  essential  to  understand 
further  the  intrinsic  dependencies  of  the  intensification 
process  on  the  vertical  eddy  diffusivity  and  on  the  surface 
drag.  In  this  article  we  focus  on  the  second  problem  and 
extend  the  study  of  Montgomery  etal  (2010)  to  include  a 
wind-speed-dependent  drag  coefficient  in  conjunction  with 
different  boundary-layer  schemes.  This  approach  permits  us 
to  assess  the  universality  of  these  dependencies.  We  use  the 
Smith  and  Thomsen  configuration  of  MM5  in  conjunction 
with  one  of  four  available  boundary-layer  schemes  in  the 
model:  the  bulk  scheme,  the  Blackadar  scheme,  the  MRF 
scheme  and  the  Gayno-Seaman  scheme. 

In  section  2  we  describe  briefly  the  numerical  model. 
The  calculations  performed  are  detailed  in  section  3 
and  supporting  interpretations  of  them  are  presented  in 
section  4.  A  discussion  of  the  results  together  with  an 
articulation  of  some  remaining  scientific  issues  is  the  subject 
of  section  5  and  a  summary  with  conclusions  is  given  in 
section  6. 

2.  The  numerical  model 

The  numerical  experiments  are  similar  to  those  described 
in  Nguyen  etal  (2008).  They  are  carried  out  using  MM5 
(version  3.6)  configured  with  three  domains:  a  coarse  mesh 
of  45  km  resolution  and  two,  two-way  nested  domains 
of  15  and  5  km  resolution,  respectively.  The  domains  are 
square  and  are  9000,  4500  and  1500  km  on  each  side.  The 
calculations  are  performed  on  an  / -plane  centred  at  20°N. 
In  all  calculations  there  are  24  a  -levels  in  the  vertical,  eleven 
of  which  are  below  1500  m  and  ten  below  1000  m  with 
diminishing  spacing  towards  the  surface.  This  number  is 
believed  to  be  adequate  for  resolving  the  boundary-layer 
dynamics  under  the  prescribed  changes  to  the  surface 
exchange  coefficients  and  vertical  mixing  of  heat  and 
momentum.  The  model  top  is  at  50  mb. 

Deep  moist  convection  is  resolved  explicitly  and 
represented  by  the  warm-rain  scheme  as  in  Montgomery 
etal  (2010).  In  addition,  we  choose  one  of  four  boundary- 
layer  schemes  available  in  the  model  as  detailed  in 
subsection  2. 1  (also  in  the  Appendix  of  Smith  and  Thomsen, 
2010).  The  warm-rain  and  boundary-layer  schemes  are 
applied  in  all  domains.  No  cumulus  parametrization  is 
used.  The  sea  surface  temperature  is  a  constant  (27  °C).  In 
contrast  to  Montgomery  et  al  (2010),  who  used  Newtonian 
relaxation  to  the  ambient  temperature  profile  as  a  radiative 
damping  scheme,  we  use  the  simple  radiative  cooling  scheme 
available  in  MM5,  which  imposes  a  temperature-dependent 
cooling  rate  on  the  order  of  l-2°Cday“h  This  choice 
helps  allay  some  concerns  that  have  been  raised  recently 
(Hakim,  2011)  on  the  appropriateness  of  the  radiative 


relaxation  scheme  that  has  been  used  in  many  idealized 
tropical  cyclone  studies  (e.g.  Rotunno  and  Emanuel,  1987; 
Persing  and  Montgomery,  2003;  Bryan  and  Rotunno,  2009). 

The  initial  vortex  is  axisymmetric  with  a  maximum 
tangential  wind  speed  of  15  m  s^^  at  the  surface  at  a  radius 
of  120  km.  The  magnitude  of  the  tangential  wind  decreases 
sinusoidally  with  height,  vanishing  at  the  top  model  level. 
The  temperature  field  is  initialised  to  be  in  gradient  wind 
balance  with  the  wind  field  using  the  method  described  by 
Smith  (2006).  The  far-field  temperature  and  humidity  are 
based  on  Jordan’s  Caribbean  sounding  for  the  hurricane 
season  (Jordan,  1958). 

For  the  purpose  of  calculating  azimuthal  averages,  the 
vortex  centre  is  defined  as  the  centroid  of  relative  vorticity 
at  900  mb  over  a  circular  region  of  200  km  radius  from  a 
‘first-guess’  centre,  which  is  determined  by  the  minimum  of 
the  total  wind  speed  at  900  mb. 

2.1.  The  boundary-layer  schemes 

The  four  boundary-layer  schemes  examined  here  are  listed 
in  T  able  1  of  Smith  and  Thomsen  (20 1 0 )  together  with  a  brief 
description  of  each.  In  the  four  main  sets  of  calculations,  the 
surface  drag  and  heat  and  moisture  exchange  coefficients 
are  modified  to  fit  the  results  of  the  Coupled  Boundary 
Layer  Air-Sea  Transfer  experiment  (CBLAST;  Black  et  al, 
2007;  Drennan  et  al,  2007;  French  et  al,  2007;  Zhang  et  al, 
2008)  to  facilitate  a  proper  comparison  of  the  schemes.  The 
surface  exchange  coefficients  for  sensible  heat  and  moisture 
are  set  to  the  same  constant,  1.2  x  10^^.  For  momentum, 
the  standard  drag  coefficient,  is  set  to 

Cdo  =  0.7x  10^^  -t-  1.4x  lO^^Jl  -  exp(-0.055|u|)}, 

where  |u|  is  the  wind  speed  at  the  lowest  model  level,  which 
is  at  about  40  m.*  This  formula  provides  an  acceptable 
approximation  to  the  data  on  Cd,  which  has  a  lot  of  scatter. 
The  drag  coefficient  and  the  exchange  coefficients  for  heat 
and  moisture  were  modified  indirectly  in  the  Blackadar, 
MRF  and  Gayno-Seaman  schemes  using  formulae  given  in 
the  Appendix  of  Hill  and  Lackman  (2009).  In  essence,  the 
friction  velocity  is  replaced  by  Vq  VCd,  where  Vq  is  the  total 
wind  speed  at  the  lowest  model  level. 

3.  The  calculations 

We  carried  out  a  series  of  numerical  experiments  of  vortex 
intensification  based  on  the  canonical  thought  experiment 
discussed  by  Nguyen  et  al  (2008),  but  using  four  of  the  five 
boundary-layer  schemes  investigated  by  Smith  and  Thomsen 
(2010).  For  all  four  schemes,  we  compare  calculations  using 
the  standard  drag  coefficient  with  those  in  which  the  drag 
coefficient  is  halved  and  doubled.  In  all  calculations,  the 
surface  exchange  coefficients  for  sensible  and  latent  heat  are 
held  constant. 


*Strictly,  the  drag  formulation  should  be  applied  at  a  height  of  10  m.  We 
have  not  adjusted  the  standard  drag  coefficient  to  the  lowest  model  level 
in  the  belief  that  the  correction  will  be  relatively  small  and  will  not  affect 
the  comparisons  between  the  calculations  described  in  section  3.  Since 
the  model  has  only  eleven  grid  points  up  to  a  height  of  about  1.5  km, 
features  on  a  vertical  scale  of  30  m  are  clearly  not  resolved.  We  have 
verified  the  adequacy  of  not  making  the  adjustment  with  one  calculation 
in  which  such  an  adjustment  was  made  and  find  that  the  differences  in 
flow  evolution  are  not  significant. 
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(a) 


(b) 


Figure  1.  Time  series  of  (a,  c,  e,  g)  maximum  azimuthally-averaged  tangential  wind,  Vniax>  and  (b,  d,  f,  h)  maximum  azimuthally-averaged  radial  inflow, 
Wmin>  iri  the  calculations  with  (a,  b)  the  bulk  scheme,  (c,  d)  the  Blackadar  scheme,  (e,  f)  the  MRF  scheme,  and  (g,  h)  the  Gayno-Seaman  scheme.  Curves  1 
(black)  denote  Cd  =  Cdo!'^^  curves  2  (blue)  Cd  =  and  curves  3  (red)  Cd  =  2Cdo- 


For  two  of  the  schemes,  we  carry  out  also  an  ensemble 
of  calculations  for  some  values  of  the  drag  coefficient  in 
which  the  initial  moisture  field  is  randomly  perturbed  with 
values  not  exceeding  0.5  gmkg^^  in  magnitude  at  altitudes 
below  the  900  mb  level  (about  1km).  The  method  follows 
that  of  Nguyen  etal.  (2008).  The  aim  is  to  demonstrate  the 
significance  of  the  differences  between  calculations  when 
the  drag  coefficient  is  halved  or  doubled.  Because  of  the 
computational  expense,  ensemble  calculations  were  carried 
out  only  for  the  Blackadar  and  Gayno-Seaman  schemes 
where  it  was  felt  that  they  were  most  necessary. 


4.  Results 

4. 1 .  Vortex  evolution 

Figure  1  shows  time  series  of  the  maximum  azimuthally- 
averaged  tangential  velocity  component,  Vmax)  and  maxi¬ 
mum  azimuthally-averaged  radial  inflow  component, 


^For  neatness  of  plotting,  we  have  defined  ‘inflow’  here  to  be  positive 
when  the  flow  is  inwards. 
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for  the  four  sets  of  calculations  described  above  including 
the  ensembles.  The  salient  points  to  note  are  as  follows. 

For  all  schemes,  both  the  maximum  tangential  wind  speed 
and  maximum  inflow  increase  steadily  with  time  as  the 
vortex  intensifies.  Both  these  quantities  show  fluctuations 
which  may  be  attributed  to  the  stochastic  nature  of  the  deep 
vortical  convective  elements  (dubbed  ‘vortical  hot  towers’  or 
VHTs)  that  are  the  fundamental  coherent  structures  driving 
the  spin-up  process  (Nguyen  et  al,  2008).  The  inflow  at  any 
given  time  is  appreciably  larger  with  the  bulk  and  Blackadar 
schemes  than  with  the  MRF  and  Gayno-Seaman  schemes. 

For  all  except  the  Gayno-Seaman  scheme,  the  inflow  is 
larger  at  most  times  when  the  value  of  the  drag  coefficient 
is  increased.  With  the  bulk  scheme,  there  are  small  time 
intervals  when,  on  account  of  the  flow  fluctuations,  Wmin  is 
almost  the  same  in  the  cases  with  Cu  —  Cpo  and  Qj  =  2Cdo. 
With  the  Blackadar  scheme,  the  ensemble  calculations 
of  Mmin  for  Cd  =  2Cdo  are  mostly  larger  than  those  for 
Cd  =  Cdo,  although  there  is  a  slight  overlap.*  With  the 
MRF  scheme,  the  differences  in  u^in  between  Cd  =  Cdo  and 
Cd  —  2Cdo  are  clear  for  much  of  the  calculation  time,  but 
they  are  not  large  and  probably  not  too  significant.  With  the 
Gayno-Seaman  scheme,  there  is  a  clear  distinction  between 
the  curves  for  Cd  —  Cdo/2  and  Cd  —  Cdo,  whereas  there  is 
a  general  overlap  of  curves  in  the  two  sets  of  ensembles  for 
Cd  —  Cdo  and  Cd  —  2Cdo>  suggesting  that  the  differences 
indicated  by  the  two  control  runs  are  not  significant.  Note 
also  that  the  spread  of  the  ensembles  for  the  Gayno-Seaman 
scheme  with  Cd  —  CdoI2  is  comparatively  small. 

The  increase  of  the  inflow  with  increasing  Cd  might  be 
anticipated,  even  if  there  were  no  change  in  the  distribution 
of  latent  heat  release  in  the  vortex  aloft.  The  reason  is  that 
increased  drag  leads  to  a  greater  reduction  of  the  near¬ 
surface  tangential  wind  speed,  at  least  in  the  outer  region 
of  the  vortex  where  there  is  subsidence  into  the  boundary 
layer.  This  greater  reduction  leads  to  a  larger  disruption  of 
gradient  wind  balance  and  a  correspondingly  larger  inward 
agradient  force  (Smith  and  Vogl,  2008;  Smith  et  al,  2009; 
Montgomery  etal,20l0)}  The  assumption  is  that  this  larger 
force  is  not  outweighed  by  an  increase  in  the  frictional  force 
in  the  radial  direction.  Of  course,  in  the  calculations  here, 
the  coupling  between  the  vortex  aloft  and  the  boundary  layer 
means  that  changes  in  the  vortex  on  account  of  changes  in 
the  pattern  of  latent  heat  release  may  explain  some  of  the 
changes  in  the  calculated  behaviour.  Such  changes  would 
be  expected  to  be  associated  mostly  with  changes  in  the 
moisture  distribution  in  the  boundary  layer  (section  5). 

In  contrast  to  the  behaviour  of  Mmim  i^max  first  increases 
as  Cd  is  increased  from  Cdo/2  to  Cdo  and  then  decreases, 
or  remains  more  or  less  the  same,  as  Cd  is  increased  further 


*One  of  the  principal  findings  of  Nguyen  etal.  (2008)  was  the 
demonstration  that  the  intensification  process  and  maximum  mean 
tangential  wind  possesses  an  intrinsic  stochastic  element  associated  with 
the  vortical  convective  elements.  In  view  of  this  result,  each  time  series 
in  Figure  1  should  be  regarded  as  a  single  member  of  an  ensemble 
generated  by,  say,  small  moisture  variations  in  the  boundary  layer. 
Based  on  the  findings  in  Nguyen  etal.  (2008),  which  show  a  spread  on 
the  order  of  10ms~*  in  the  maximum  azimuthaUy-averaged  tangential 
velocity  amongst  ensemble  members  during  the  intensification  phase, 
the  differences  between  runs  that  are  less  than  this  amount  would  not  be 
regarded  as  significant. 

^The  agradient  force  is  the  sum  of  the  radial  pressure  gradient  force  and 
the  Coriolis  and  centrifugal  forces,  usually  expressed  per  unit  mass  of 
air. 


to  2Cdo-  In  part,  the  reason  is  that  the  increased  boundary- 
layer  inflow  brings  air  parcels  closer  to  the  vortex  centre 
leading  to  the  tendency  for  them  to  spin  faster  due  to  the 
partial  conservation  of  absolute  angular  momentum.^  This 
tendency  is  offset  by  the  frictional  torque,  which  reduces  the 
absolute  angular  momentum.  If,  on  account  of  an  increase 
in  Cd,  air  parcels  converge  sufficiently  rapidly,  the  loss  of 
absolute  angular  momentum  may  be  more  than  offset  by  the 
radial  displacement  towards  the  centre,  leading  to  a  higher 
tangential  wind  speed  maximum  (Smith  and  Vogl,  2008; 
Smith  et  al,  2009;  Montgomery  et  al,  2010).  Figure  1  shows 
that  the  circumstances  when  the  convergence  of  absolute 
angular  momentum  dominates  is  in  the  lower  range  of  drag 
coefficients  studied:  between  Cdo/2  and  Cdo-  However,  as 
Cd  is  increased  further,  the  increased  torque  wins.  Thus, 
all  schemes  have  an  optimum  drag  for  the  maximum 
intensification  of  Vmax.  This  result  is  to  be  expected  from 
the  results  of  Montgomery  etal  (2010),  since  the  case  of 
zero  drag  gives  no  sustained  intensification  and  the  case  of 
extreme  drag  gives  little  intensification  also. 

The  foregoing  results  appear  to  he  significant  in  the 
light  of  Emanuel’s  theory  for  the  potential  intensity  of 
a  tropical  cyclone  (Emanuel,  1986;  Emanuel,  1995;  Bister 
and  Emanuel,  1998),  which  predicts  a  decrease  in  v^ax 
as  the  drag  coefficient  increases  (e.g.  Emanuel,  1995,  his 
Figure  1).  The  three-dimensional  calculations  here  and  in 
Montgomery  etal  (2010)  do  not  exhibit  such  behaviour 
as  the  drag  coefficient  is  varied.  Contrary  to  a  prediction 
of  Emanuel’s  (1995)  theory,  there  is  an  optimum  value  of 
Q/Cd  for  the  mature  intensity  if  the  enthalpy  exchange 
coefficient  is  held  fixed.  At  this  stage  we  are  unable  to 
provide  a  satisfactory  explanation  for  this  difference  in 
behaviour.  Although  Emanuel’s  theory  is  for  the  purported 
steady  state  and  our  calculations  are  terminated  after  five 
days,  our  calculations  do  appear  to  be  approaching  at  least 
a  quasi-steady  state.  One  possibility  for  the  discrepancy 
might  be  the  oversimplified  representation  of  the  boundary 
layer  in  Emanuel’s  theory,  or  his  use  of  an  axisymmetric 
model.  Although  the  calculations  of  Bryan  and  Rotunno 
(2009)  largely  support  Emanuel’s  theoretical  predictions 
(their  Figure  6),  they  used  an  axisymmetric  model  also. 
Results  of  a  very  recent  paper  by  Bryan  (2012)  support  the 
findings  of  Montgomery  et  al  (2010)  and  those  herein  for  the 
intensification  stage,  but  agree  with  Emanuel’s  predictions 
for  the  maximum  intensity.  However,  it  is  unclear  whether 
Bryan’s  calculations  have  reached  a  steady  state  also.  In  fact, 
as  discussed  in  section  5.2,  it  is  unclear  whether  a  strict 
steady  state  can  be  realized  in  which  all  time  tendencies  are 
identically  zero  globally.  We  believe  that  further  research  on 
this  issue  is  called  for. 

4.2.  Low-level  wind  structure 

Figure  2  shows  radius-height  cross-sections  of  the 
azimuthally-averaged  radial  velocity  component,  u,  and 
tangential  velocity  component,  v,  in  the  lowest  3  km, 
averaged  during  the  period  1 14-120  h,  for  the  calculations 


^The  absolute  angular  momentum,  M,  is  defined  by  rv  -F  j/r^,  where  r 
is  the  radius,  v  is  the  tangential  wind  component  and  /  is  the  Coriolis 
parameter.  Near  the  surface,  M  is  not  conserved  on  account  of  the 
frictional  torque.  Nevertheless,  the  convergence  of  M  leads  to  spin-up  if 
air  parcels  can  be  brought  to  small  radii  quickly  enough,  before  losing 
an  appreciable  amount  of  M  (Smith  and  Vogl,  2008;  Smith  et  al.,  2009). 
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Figure  2.  Radius-height  cross-sections  of  azimuthally-averaged  radial  (thin  lines  with  negative  values  dashed)  and  tangential  (bold  lines)  wind  speed 
components  in  the  lowest  3  km  averaged  over  the  period  1 14-120  h  for  (a,  c,  e)  the  Blackadar  scheme  and  (b,  d,  f)  MRF  scheme.  (a,z,b)  show  Co  =  Cdi)/2, 
(c,  d)  Co  =  Coo,  and  (e,  f)  Co  =  2Cdo-  The  contour  interval  is  5ms“'  for  both  components,  except  for  the  radial  component  in  (b),  where  it  is  only 
1  m  s“* .  This  figure  is  available  in  colour  online  at  wileyonlinelibrary.com/journal/qj 


with  the  Blackadar  and  MRF  schemes.  The  three  panels  for 
each  scheme  show  fields  for  Co  =  CdoI^,  Cdo,  and  2Cdo- 
Prominent  features  are  that  Vmax  occurs  at  low  levels  near 
the  top  of  the  frictionally  induced  inflow  layer  in  aU  cases, 
while  Mmin  occurs  very  close  to  the  surface,  in  accord  with  the 
findings  of  Smith  et  al.  (2009).  However  there  are  significant 
differences  in  the  behaviour  of  the  two  schemes  to  increasing 
Cd  as  follows. 

4.2. 1 .  Blackadar  scheme 

With  the  Blackadar  scheme,  the  radius  of  Vmax  decreases 
significantly  as  Co  increases,  being  about  80  km  for  Co  — 
CoojT-^  40  km  for  Co  —  Coo,  and  33  km  for  Co  —  2Coo- The 
radius  of  Wmin  decreases  also  as  Co  is  increased.  A  physical 
explanation  for  this  behaviour  follows  directly  from  the  ideas 
discussed  in  the  previous  section.  The  increased  friction  leads 
to  a  larger  gradient  wind  imbalance  in  the  boundary  layer. 
With  the  higher  Co  value,  inflowing  rings  of  boundary-layer 
air  are  subjected  to  a  greater  agradient  force  and  hence  are 
driven  farther  inwards  before  rising  out  of  the  boundary 
layer  and  ascending  into  the  eyewall  updraught.  The  end 
result  is  an  enhanced  maximum  tangential  wind,  despite 
some  additional  loss  of  absolute  angular  momentum  en 


route  (Nguyen  et  al,  2002;  Smith  et  al,  2009;  Montgomery 
et  al,  2010). 

The  foregoing  argument  assumes  that  the  boundary-layer 
depth,  as  characterized  by  the  depth  of  appreciable  inflow,  I' 
remains  essentially  unaltered  as  the  drag  coefficient  changes. 
However,  as  noted  by  Smith  and  Vogl  (2008),  the  effects 
of  surface  drag  are  distributed  through  the  depth  of  the 
boundary  layer  and  a  change  in  the  depth  through  which  the 
net  radial  force  acts  will  change  the  strength  of  the  inflow. 
As  seen  in  Figure  2,  the  depth  of  the  inflow  increases  with 
increasing  drag,  but  only  marginally  as  Co  goes  from  Coo 
to  2Coo-  This  behaviour  is  consistent  with  the  predictions 
of  linear  theory,  which  shows  that  the  height  of  maximum 


II  Following  Smith  et  al.  (2009),  we  use  the  term  boundary  layer  to  describe 
the  shallow  layer  of  strong  inflow  near  the  sea  surface  that  is  typically 
500  m  to  1  km  deep  and  which  arises  largely  because  of  the  frictional 
disruption  of  gradient  wind  balance  near  the  surface.  Although  there 
is  some  inflow  throughout  the  lower  troposphere  associated  with  the 
balanced  response  of  the  vortex  to  latent  heat  release  in  the  eyewall  clouds 
(Bui  et  al,  2009),  the  largest  radial  wind  speeds  are  confined  within  the 
lowest  kilometre  and  delineate  clearly  the  layer  in  which  friction  effects 
are  important  (i.e.  where  there  is  gradient  wind  imbalance;  Figure  6  of 
Smith  et  al,  2009)  from  the  region  above  where  they  are  not. 
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inflow  generally  increases  with  increasing  Co  (Figure  1  of 
Kepert,  2001). 

Vogl  and  Smith  (2009)  presented  a  scale  analysis  of  the 
boundary  layer  showing  that  the  boundary-layer  depth  is 
proportional  to  the  square  root  of  the  eddy  diffusivity.  While 
the  eddy  diffusivity  is  assumed  constant  in  that  analysis,  in 
more  general  situations  the  boundary-layer  depth  may  be 
assumed  to  increase  with  some  vertical  average  of  the  eddy 
diffusivity.  The  distribution  of  eddy  diffusivity  predicted  by 
the  various  schemes  is  discussed  further  in  section  4.3. 

4.2.2.  MRF  scheme 

With  the  MRF  scheme,  the  layer  of  appreciable  inflow  is 
much  deeper  than  with  the  Blackadar  scheme,  a  result  that 
is  consistent  with  the  relatively  large  eddy  diffusivity  in  this 
scheme.  We  show  in  section  4.3  that  the  eddy  diffusivities 
predicted  by  the  MRF  scheme  are  much  larger  than  for 
the  Blackadar  scheme  for  the  same  representation  of  drag 
coefficient.  There  is  again  a  relatively  large  change  in  the 
radius  ofvmax  as  Cd  goes  from  CdoI'2-  to  Cdo,  butlittle  change 
as  it  goes  from  Cdo  to  2Cdo-  Moreover,  as  for  the  Blackadar 
scheme,  the  radius  of  Wmin  decreases  as  Cd  goes  from  CdoI2 
to  Cdo  but,  unlike  the  behaviour  with  the  Blackadar  scheme, 
it  increases  as  Cd  goes  from  Cdo  to  2Cdo-  It  follows  that,  for  an 
increase  in  drag  coefficient  beyond  the  standard  value  in  the 
MRF  scheme,  the  tangential  wind  tendency  accompanying 
the  increased  convergence  of  absolute  angular  momentum 
is  much  more  than  compensated  by  the  negative  tendency 
of  friction.  The  weaker  inflow  with  the  MRF  scheme  has 
interesting  consequences  for  the  pattern  of  vertical  motion 
as  discussed  in  section  5.4. 

4.2.3.  Other  schemes 

For  space  reasons  we  have  not  shown  figures  corresponding 
to  Figure  2  for  the  other  two  schemes,  but  simply  summarize 
the  results.  The  changes  in  low-level  structure  with  the 
change  in  drag  coefficient  in  the  case  of  the  bulk  scheme  is 
similar  to  that  of  the  Blackadar  scheme,  but  the  vortex  in  the 
case  with  Cd  =  CdoI2  is  more  intense  than  its  counterpart 
with  the  Blackadar  scheme  (compare  Figure  2(a)  and  (c)). 
The  behaviour  in  the  case  of  the  Gayno-Seaman  scheme  is 
similar  to  that  with  the  Blackadar  scheme  in  that  the  vortex 
in  the  1 14-120  h  average  with  Cd  =  CdoI2  is  considerably 
weaker  than  those  for  Cd  and  Cd  —  2Cdo-  However,  the 
vortex  for  Cd  —  2Cdo  is  only  slightly  weaker  than  that  for 
Cd  —  Cdo- 

4.3.  Eddy  diffusivity 

Some  of  the  differences  in  behaviour  between  schemes  can 
be  interpreted  in  terms  of  the  distribution  of  vertical  eddy 
diffusivity  that  they  produce.  While  the  schemes  investigated 
have  various  degrees  of  sophistication,  almost  all  seek  to 
determine  some  local  value  of  turbulent  diffusivity,  K{r,  z), 
to  close  the  momentum  and  thermodynamic  equations.  An 
exception  is  the  Blackadar  scheme  in  the  (dry)  convective 
regime,  which  uses  a  non-local  mixing  algorithm.  In  our 
calculations,  we  found  that  this  algorithm  is  not  invoked 
within  200  km  of  the  centre  of  the  vortex  circulation.  In 
other  stability  regimes  of  this  scheme,  the  formulation  for 
K{r,z)  is  exactly  the  same  as  in  the  bulk  scheme.  In  both 
these  schemes,  K  is  the  product  of  the  square  of  a  mixing 


length,  the  vertical  wind  shear,  and  a  specified  function  of  the 
Richardson  number  (Blackadar,  1976).  In  the  MRF  scheme, 
the  determination  ofi<C  is  based  on  empirical  formulae,  while 
in  the  Gayno-Seaman  scheme  it  is  related  to  the  turbulent 
kinetic  energy,  which  is  a  prognostic  quantity  in  this  scheme. 

Figure  3  shows  radius-height  cross-sections  of  the  vertical 
eddy  diffusivity  in  the  lowest  3  km,  averaged  during  the 
period  1 14-120  h,  for  the  calculations  with  the  Blackadar 
and  MRF  schemes  for  Cd  —  CdoI2,  Cdo  and  2Cdo-  The 
patterns  in  all  four  schemes  including  the  two  not  shown 
are  similar,  and  similar  also  to  those  found  by  Braun 
and  Tao  (2000),  Smith  and  Thomsen  (2010)  and  Kepert 
(2012).  In  particular,  there  is  a  low-level  maximum  of  K 
in  the  region  of  maximum  tangential  wind  speed  as  well 
as  elevated  values  in  the  eyewall  updraught.  The  main 
updraught  is  indicated  in  the  figure  by  the  0.2  m  s^'  contour 
of  vertical  velocity  (0.1  ms^^  contour  in  (b)).  There  is  a 
strong  radial  gradient  of  K  also  and,  for  the  two  schemes 
in  the  figure,  the  values  of  K  increase  with  increasing  drag 
coefficient. 

Maximum  values  of  K  for  all  the  schemes  are  shown  as 
time  series  in  Figure  4,  where  it  is  seen  that  the  MRF  scheme 
is  by  far  the  most  diffusive  with  a  maximum  of  nearly 
600  m^  s^'  (note  the  five-fold  increase  in  contour  interval 
needed  for  this  scheme).  This  result  explains  why  the  inflow 
layer  is  much  deeper  in  the  MRF  scheme  and,  because  the 
surface  drag  is  therefore  distributed  over  a  much  deeper 
layer,  why  the  local  disruption  of  gradient  balance  is  smaller 
leading  to  a  substantially  weaker  inflow. 

Maximum  values  of  K  with  the  bulk  scheme  are 
comparable  with,  but  systematically  a  little  larger  than, 
those  with  the  Blackadar  scheme  and  lie  in  the  range 
60-120  m^s“k  In  contrast,  the  Gayno-Seaman  scheme  is 
relatively  diffusive  with  maximum  values  of  K  reaching 
about  250m^s~^  All  schemes  except  the  Gayno-Seaman 
scheme  predict  a  marked  increase  in  the  maximum  K  as 
the  drag  coefficient  is  increased.  With  the  Gayno-Seaman 
scheme,  the  maximum  K  changes  relatively  little,  although 
there  are  noticeable  differences  in  the  distribution  of  K  when 
the  drag  coefficient  is  increased  (not  shown).  For  example, 
the  areal  extent  within  the  K=50m^s^^  contour  in 
the  azimuthal-mean,  radius-height  cross-section  increases 
significantly,  but  not  as  much  as  with  the  MRF  scheme. 

The  only  observational  estimates  for  K  that  we  are 
aware  of  are  those  analysed  recently  from  flight-level  wind 
measurements  at  an  altitude  of  about  500  m  in  hurricanes 
Allen  {1980}  and  Hugo  (1989)  by  Zhang  et  al.  {2011). In  Hugo 
these  were  about  llOm^s^^  beneath  the  eyewall,  where 
the  near-surface  wind  speeds  were  about  60ms^\  and  in 
Allen  they  were  up  to  74  m^  s^\  where  wind  speeds  were 
about  72  m  s^^  (Zhang  et  al,  2011).  One  might  be  tempted 
to  judge  that  the  MRF  and  Gayno-Seaman  schemes  are 
much  too  diffusive,  whereas  the  other  schemes  have  broadly 
realistic  diffusivities,  but  it  would  be  premature  to  draw  firm 
conclusions  from  a  comparison  with  only  two  observational 
estimates! 

5.  Discussion  and  some  remaining  scientific  issues 

The  interpretations  given  above  for  the  vortex  behaviour 
as  a  result  of  changing  the  boundary-layer  scheme,  or 
changing  the  drag  coefficient  within  a  scheme,  go  some  way 
to  providing  basic  understanding  of  the  problem.  However, 
they  fall  short  of  providing  a  complete  theory,  which  would 
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Figure  3.  Radius-height  cross-section  of  the  time-mean,  azimuthally-averaged  eddy  diffusivity,  K{r,z),  during  the  period  114-120h  for  (a,c,e)  the 
Blackadar  scheme  (with  contour  interval  10m^s“*)  and  (b,d,f)  MRF  scheme  (with  contour  interval  50m^s“*).  (a,b)  Cd  =  standard/2;  (c.d)  Cd  = 
standard;  (e,  f)  Cd  =  2x  standard.  The  main  updraught  is  indicated  in  the  figure  by  the  bold  black  contour,  which  depicts  the  0.2  m  s“'  contour  of  the 
same  time-mean  azimuthally-averaged  vertical  velocity,  except  in  (b)  where  the  contour  value  is  0.1  ms“*.  This  figure  is  available  in  colour  online  at 
wileyonlinelibrary.com/journal/qj 


Figure  4.  Time  series  of  maximum  azimuthally-averaged  eddy  diffusivity,  ffniax>  for  the  four  modified  boundary-layer  schemes:  (a)  bulk  (dashed)  and 
Blackadar  (solid)  schemes;  and  (b)  MRF  (dashed)  and  Gayno-Seaman  (solid)  schemes.  Curves  1  refer  to  Cd  =  Cdo/2,  curves  2  to  Cd  =  Cdo,  and  curves 
3  to  Cd  =  Cdo.  Note  the  different  scales  along  the  ordinate.  This  figure  is  available  in  colour  online  at  wileyonlinelibrary.com/journal/qj 


entail  the  construction  of  a  physically  consistent  model 
to  mimic  the  behaviour  shown.  So  far,  the  interpretations 
have  focussed  on  boundary-layer  processes,  but  a  complete 
theory  would  require  consideration  of  processes  above  the 
boundary  layer  as  well.  Even  though  the  enthalpy  exchange 


coefficients  are  held  constant  in  all  calculations,  a  change  in 
the  surface  wind  speed  brought  about  by  a  change  in  the  value 
of  Cd  will  lead  to  a  change  in  the  surface  enthalpy  flux,  unless 
the  air-sea  disequilibrium  were  to  exactly  compensate.  The 
latter  scenario  is  highly  improbable.  A  change  in  enthalpy 
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fluxes  (primarily  latent  heat)  must  feed  back  to  affect  the 
spatial  distribution  of  diabatic  heating  and  the  aggregate 
effects  of  buoyancy  associated  with  deep  convection.  At  this 
stage  we  are  unaware  of  a  satisfactory  theory  to  link  these 
changes  with  a  change  in  Cd-  The  most  we  can  do  here  is  to 
attempt  to  articulate  the  main  challenges  in  developing  such 
a  theory. 

First  we  note  that  the  interpretations  given  herein 
are  based  on  axisymmetric  reasoning.  While  the  flow 
asymmetries  in  the  form  of  rotating  deep  convective  clouds** 
are  quantitatively  important  in  the  flow  dynamics,  it  is  their 
collective  effect  which  is  relevant  to  the  present  discussion, 
so  that  axisymmetric  reasoning  should  suffice. 

5.1.  The  intensification  stage 

We  consider  first  the  intensification  stage.  As  shown  in  Bui 
etal.  (2009),  the  convective  clouds  provide  an  aggregate 
heating  rate  that  drives  a  system-scale  convergence  of 
absolute  angular  momentum  above  the  boundary  layer, 
where  this  quantity  is  approximately  materially  conserved. 
The  net  effect  of  this  convergence  is  to  spin  up  the  circulation 
there.  The  situation  in  the  boundary  layer  is  more  subtle.  We 
would  argue  that  the  vortex-scale  inflow  in  the  boundary 
layer  is  a  response  to  the  azimuthally-averaged  tangential 
wind  at  the  top  of  the  boundary  layer  and,  because  of 
the  parabolic  nature  of  the  boundary-layer  equations,^^ 
convection  cannot  ‘suck’  inflow  in  the  boundary  layer, 
except  possibly  near  to  convection  where  boundary-layer 
theory  breaks  down.  As  noted  earlier,  the  inflow  in  the 
boundary  layer  is  much  stronger  than  that  above  on  account 
of  the  inwards-directed  agradient  force,  but  absolute  angular 
momentum  in  this  layer  is  lost  by  the  frictional  torque. 
Typically,  in  the  inner-core  region,  the  tangential  wind 
becomes  supergradient  and,  as  a  result,  the  maximum 
tangential  wind  speed  occurs  within,  but  near  the  top  of,  the 
boundary  layer. 

The  diabatic  heating  rate  associated  with  deep  convec¬ 
tion  depends  on  both  the  vertical  velocity  and  specific 
entropy  of  ascending  air,  the  radial  distribution  of  which 
must  he  determined  by  the  boundary-layer  dynamics  and 
thermodynamics.  The  boundary-layer  dynamics  must  deter¬ 
mine  also  the  radial  distribution  of  horizontal  momentum 
(including  angular  momentum)  of  the  ascending  air  follow¬ 
ing  its  spin-up  in  the  boundary  layer.  Above  the  boundary 
layer,  a  scale  analysis  indicates  that  the  flow  is  close  to 
gradient  wind  balance  (Willoughby,  1979),  so  that  the  sec¬ 
ondary  circulation  is  governed  there  by  the  Sawyer-Eliassen 
equation  (Shapiro  and  Willoughby,  1982;  Bui  et  al,  2009), 
at  least  during  periods  of  intensity  change  (see  below).  In 
the  absence  of  friction,  the  forcing  term  in  this  equation  is 
primarily  the  radial  gradient  of  diabatic  heating.  However, 
there  is  an  additional  forcing  at  the  lower  boundary  (i.e.  the 
top  of  the  frictional  boundary  layer)  on  account  of  air  being 
expelled  from  or  drawn  into  the  boundary  layer. 


**An  examination  of  the  vertical  motion  fields  at  the  850  mb  level  shows 
that  the  character  of  the  convection  does  not  change  in  any  obvious  way 
with  a  change  in  boundary-layer  scheme. 

^^The  derivation  of  the  boundary-layer  equations  for  a  tropical  cyclone 
is  given  in  Smith  (1968)  and  Vogl  and  Smith  (2009)  on  the  basis  of  a 
scale  analysis  of  the  Navier-Stokes  equations. 


5.2.  The  mature  stage 

We  consider  now  the  mature  stage  of  development.  In  a  strict 
steady  state,  the  steps  taken  to  derive  the  Sawyer-Eliassen 
equation  become  degenerate  because  the  time  tendencies 
that  are  eliminated  in  the  derivation  are  identically  zero. 
Then,  in  the  absence  of  friction  and/or  eddy  fluxes  of 
tangential  momentum  and  moist  entropy,  the  azimuthally- 
averaged  secondary  circulation  above  the  boundary  layer 
must  be  along  both  absolute  angular  momentum  surfaces 
and  moist  isentropic  surfaces,  i.e.  these  two  sets  of  surfaces 
must  be  congruent.  However,  these  constraints  are  not 
sufficient  by  themselves  to  determine  the  flow,  even  if  the 
flow  is  in  balance.^^  Indeed,  as  discussed  in  the  footnote 
here,  it  is  unclear  that  a  steady  solution  of  even  the  balance 
problem  exists  globally.  This  issue  presents  a  vexing  problem 
in  relating  the  change  in  diabatic  heating  distribution  with  a 
change  in  the  drag  coefficient  in  the  mature  stage. 

Two  further  problems  of  constructing  a  steady-state 
theory  for  a  hurricane  are  the  fact  that  the  inner-core 
boundary  layer  is  not  in  gradient-wind  balance  and  the  top 
of  the  boundary  layer  is  open.  At  radii  where  there  is  ascent 
out  of  the  boundary  layer,  an  open  boundary  condition 
is  required  to  allow,  inter  alia,  the  high  tangential  wind 
speeds  that  are  spun  up  in  the  boundary  layer  to  be  lofted 
into  the  interior  flow.  The  formulation  of  such  a  condition 
is  problematic  as  discussed  by  Smith  and  Montgomery 
(2010)  and  it  seems  wrong  to  us  to  constrain  the  flow 
to  return  to  a  prescribed  gradient  wind  as  is  typically 
assumed  in  steady  boundary-layer  theory.  Because  of  the 
tight  coupling  between  the  boundary  layer  and  the  interior 
flow,  we  would  be  uncomfortable  about  investigating  the 
effects  of  changing  the  drag  coefficient  by  focussing  only 
on  the  steady  boundary-layer  response  and  believe  that  it  is 
essential  to  solve  for  the  vortex  flow  as  a  whole. 

5.3.  Effects  of  changing  Cd  on  the  primary  circulation 

Notwithstanding  the  foregoing  issues,  it  may  be  of  interest 
to  the  reader  to  see  the  differences  in  the  flow  above  the 
boundary  layer  as  a  result  of  changes  in  the  drag  coefficient. 
For  this  reason,  we  show  in  Figures  5  and  6  some  details  of 
the  flow  above  the  boundary  layer  for  the  four  boundary- 
layer  schemes.  We  show  these  figures  to  provide  a  context 
for  future  theoretical  investigations,  but  at  this  time  we  are 
unable  to  provide  a  satisfactory  interpretation  for  the  range 
of  behaviour  surveyed.  Figure  5  shows  radial  profiles  of 
azimuthally-averaged  tangential  wind  speed,  Vg,  at  a  height 
of  2  km  for  the  four  schemes  considered  and  for  the  three 


more  strict  constraint  on  steady  balance  fiow  may  be  obtained  from 
the  geopotential  tendency  equation  derived  by  Shapiro  and  Montgomery 
(1993)  and  Vigh  and  Schubert  (2009).  If  the  flow  is  symmetrically  stable 
this  second-order  linear  partial  differential  equation  is  elliptic  with  a 
forcing  term  that  is  proportional  to  the  product  of  potential  vorticity, 
P,  with  the  6  derivative  of  6 /P  along  an  absolute  momentum  surface, 
where  9  is  the  potential  temperature  and  0  is  its  material  derivative.  Vigh 
and  Schubert  refer  to  this  product  as  the  cyclogenesis  function.  Unlike  the 
Sawyer-Eliassen  equation,  the  derivation  of  the  geopotential  tendency 
equation  is  not  degenerate  for  steady  flow,  implying  that,  for  such  a 
flow,  the  cyclogenesis  function  must  be  identically  zero.  This  condition 
determines  the  secondary  circulation  in  terms  of  the  primary  circulation 
and  the  diabatic  heating  (and,  in  general,  other  forcing  terms),  and 
although  the  requirement  of  continuity  imposed  a  further  constraint,  it 
does  not  allow  for  solution  for  the  primary  circulation  and  mass  field. 
Indeed,  it  is  even  unclear  that  a  global  steady-state  solution  exists. 
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Figure  5.  Radial  profiles  of  time-mean  azknuthally-averaged  tangential  wind  speed  at  a  height  of  2  km  for  (a)  the  bulk  and  Blackadar  schemes,  and  (b) 
the  MRF  and  Gayno-Seaman  schemes  during  the  period  1 14-120 h.  Curves  1, 2, 3  refer  to  runs  with  Co  =  standard/2;  Co  =  standard;  and  Co  =  2x 
standard,  respectively.  This  figure  is  available  in  colour  online  at  wileyonlinelibrary.com/journal/qj 


Figure  6.  Radius-height  cross-sections  of  azimuthally-averaged  vertical  velocity  for  (a,  c,  e)  the  Blackadar  scheme,  and  (b,  d,  f)  the  MRF  scheme  time 
averaged  between  114  and  120  h,  for  (a,  b)  Co  =  standard/2,  (c,  d)  Co  =  standard,  and  (e,  f)  Cd  =  2x  standard.  Thin  contours  are  at  0,  0.1,  0.2,  0.3  and 
0.4  m  s“*  and  bold  contours  are  at  0.5  m  s“*  intervals  thereafter.  This  figure  is  available  in  colour  online  at  wileyonlinelibrary.com/journal/qj 


values  of  drag  coefficient  for  each  scheme.  The  wind  speed 
Vg  should  be  a  reasonable  approximation  to  the  gradient 
wind  at  the  top  of  the  boundary  layer,  but  not  to  the  total 
surface  wind,  which  is  one  factor  determining  the  latent  heat 
flux. 

A  common  feature  of  all  the  schemes  is  that,  like  Vmax) 
the  maximum  value  of  Vg,  is  largest  for  Cp  —  Cdo- 


However,  unlike  v^ax  (cf.  Figure  1),  for  the  bulk,  Blackadar 
and  Gayno-Seaman  schemes,  is  larger  for  Co  =  CdoI2 
than  for  Cd  =  2Cdo-  The  exception  is  the  MRF  scheme, 
where  Vg^^  is  smallest  for  Cd  —  Cdo/2.  The  radius  of 
maximum  Vg,  say  rvgm^xt  varies  also  with  Cd-  With  the 
Blackadar  and  MRF  schemes,  r^gmax  is  about  twice  as  large 
for  Cd  —  CdoI2  as  with  Cd  —  Cdo^  but  the  difference  is  less 
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for  the  bulk  scheme  and  much  less  for  the  Gayno-Seaman 
scheme.  With  all  except  the  bulk  scheme,  rvgm^  decreases 
monotonically  with  increasing  Cd,  but  the  radial  profile 
of  Vg  with  the  bulk  scheme  with  Cd  =  2Cdo  is  rather 
flat. 

5.4.  Effects  of  changing  Cd  on  the  secondary  circulation 

A  feel  for  the  differences  in  the  secondary  circulation  as  Cd 
is  halved  or  doubled  is  provided  by  vertical  cross-sections  of 
the  azimuthaUy-averaged  vertical  velocity  for  the  Blackadar 
and  MRF  schemes  shown  in  Figure  6.  For  both  schemes,  the 
mean  updraught  that  marks  the  time- averaged  ascent  in  the 
eyewaU  increases  significantly  in  strength  as  Cd  increases 
from  CdoI2  to  Cdo-  Also,  the  radius  of  the  maximum 
updraught  contracts  inwards.  For  the  Blackadar  scheme,  the 
updraught  strength  increases  further  as  Cd  is  increased  to 
2Cdo  and  there  is  a  further  contraction  in  radius.  A  striking 
feature  of  this  scheme  for  the  two  larger  values  of  Cd  is 
the  pronounced  low-level  maximum  that  is  a  consequence 
of  the  rapid  deceleration  of  the  boundary-layer  inflow  and 
its  eruption  into  the  updraught.  This  eruption  corresponds 
with  a  separation  of  the  boundary  layer  brought  about  by 
inertial  (i.e.  Coriolis  and  centrifugal)  forces  in  conjunction 
with  vertical  perturbation  pressure  gradient  forces^^  as 
the  tangential  wind  component  becomes  supergradient 
(e.g.  Smith  et  al,  2009;  Smith  and  Montgomery,  2010). 
Subsequent  to  this  eruption,  the  updraught  decelerates  on 
account  of  the  static  stability  until  cloud  buoyancy  becomes 
sufficient  to  re-accelerate  it.  For  the  MRF  scheme,  the  vertical 
velocity  field  changes  more  in  pattern  than  in  strength  as 
Cd  is  increased  to  2Cdo  and  there  is  little  radial  contraction. 
Notably  the  maximum  vertical  velocity  is  lowered  and  there 
is  no  low-level  maximum  as  in  the  Blackadar  scheme.  We 
attribute  the  absence  of  the  low-level  maximum  with  the 
MRF  scheme  to  the  much  higher  vertical  diffusivity  of  the 
scheme,  which  tends  to  mitigate  the  inertial  dynamics  of  the 
boundary  layer. 

5.5.  Effects  of  changing  Cd  on  the  surface  fluxes 

It  has  been  long  recognized  that  the  latent  heat  supplied  by 
the  ocean  is  the  primary  energy  source  for  a  tropical-cyclone 
(e.g.  Riehl,  1954).  For  this  reason,  theories  for  tropical- 
cyclone  behaviour  have  laid  great  emphasis  on  incorporating 
these  fluxes  in  the  physical  formulation.  A  widely  accepted 
paradigm  for  tropical-cyclone  intensification  is  based  on 
the  supposed  feedback  between  the  surface  wind  speed  and 
the  latent  heat  flux  (Rotunno  and  Emanuel,  1987).  While 
this  feedback  process  has  been  called  into  question  recently, 
it  remains  true  that  modest  sea-to-air  vapour  fluxes  are 
essential  for  supplying  buoyancy  to  rotating  deep  convection 
(Montgomery  et  al,  2009).  With  this  background  we  show 
here  the  effect  of  changing  the  drag  coefficient  on  the 
radial  profiles  of  azimuthally- averaged  latent  heat  flux  in  the 
mature  stage  of  the  calculation. 

Figure  7  shows  radial  profiles  of  azimuthally-averaged 
latent  heat  flux  at  the  ocean  surface  in  the  mature  stage. 


^  ^  Vertical  perturbation  pressure  gradient  forces  are  required  to  accelerate 
air  parcels  upwards  as  they  are  slowed  down  in  the  radial  direction. 
As  these  forces  become  appreciable,  the  classical  boundary-layer 
approximation  that  the  transverse  pressure  gradient  in  the  boundary 
layer  is  equal  to  its  free  stream  value  is  no  longer  valid. 


These  profiles  exhibit  a  similar  behaviour  for  all  schemes 
as  Cd  is  increased,  with  the  maximum  flux  occurring 
when  Cd  —  Cdo  and  the  minimum  when  Cd  —  CdoI2.  This 
behaviour  mirrors  that  of  Vmaxj  but  not  that  of  V^,  because 
the  latter,  in  particular,  does  not  characterize  the  profile  of 
total  surface  wind,  and  even  if  it  did,  the  flux  depends  also  on 
the  thermodynamic  disequilibrium  across  the  sea  surface. 
The  radius  of  the  maximum  flux  generally  decreases  with 
increasing  Cd,  which  is  presumably,  in  part,  a  reflection 
of  the  contraction  of  the  maximum  in  both  wind  speed 
components  near  the  surface  (not  shown). 

At  this  stage  it  does  not  seem  possible  to  give  a  complete 
interpretation  of  the  changes  in  the  latent  heat  flux  with  the 
changes  in  Cd,  although  the  results  suggest  that  the  fluxes 
are  largely  slaved  to  the  boundary-layer  dynamics  through 
the  distribution  of  surface  wind  speed.  A  missing  link  in 
providing  a  complete  theory  is  knowing  how  to  relate  the 
surface  fluxes  to  local  cloud  buoyancy,  which  is  necessary  in 
an  aggregate  sense  to  drive  a  deep  overturning  circulation. 

6.  Summary  and  conclusions 

We  have  carried  out  idealized  three-dimensional  numerical 
simulations  to  investigate  the  sensitivity  of  tropical-cyclone 
intensification  to  changes  in  the  surface  drag  coefficient 
in  the  prototype  intensification  problem  discussed  by 
Nguyen  et  al  (2008).  Changing  the  drag  coefficient  provides 
insight  into  unbalanced  effects  in  the  boundary  layer  and 
their  impact  on  the  vortex  evolution.  It  provides  also 
further  understanding  of  the  intrinsic  dependencies  of  the 
intensification  process  on  the  vertical  eddy  diffusivity  and 
on  the  surface  drag.  Because  of  the  tight  coupling  between 
the  boundary  layer  and  the  interior  flow,  we  have  purposely 
solved  the  vortex  problem  as  a  whole  and  not  just  the 
boundary-layer  response  to  the  drag  coefficient.  The  study, 
which  uses  the  MM5  model,  extends  that  of  Montgomery 
etal  (2010)  to  include  a  wind-speed-dependent  drag 
coefficient  in  conjunction  with  four  different  boundary-layer 
schemes  available  in  MM5:  the  bulk  scheme,  the  Blackadar 
scheme,  the  MRF  scheme  and  the  Gayno-Seaman  scheme. 
To  facilitate  direct  comparison,  the  schemes  were  modified 
slightly  to  have  the  same  drag  coefficient  formulation  and 
the  same  constant  exchange  coefficients  for  sensible  heat 
and  moisture. 

For  a  standard  value  of  the  drag  coefficient  consistent 
with  recent  observations,  we  have  investigated  the  changes 
in  vortex  evolution  and  mature  intensity  when  this  standard 
value  is  halved  or  doubled.  The  intensity  is  characterized  by 
the  maximum  azimuthally-averaged  tangential  wind  speed. 
We  find  that  the  change  in  intensity  when  the  drag  coefficient 
is  halved  or  doubled  from  its  standard  value  depends  on 
the  boundary-layer  scheme  employed.  In  all  schemes,  the 
intensity  increases  as  the  drag  coefficient  is  increased  from 
half  its  standard  value  to  the  standard  value  and  then 
declines  as  the  drag  coefficient  is  doubled  from  its  standard 
value.  Thus  the  standard  drag  coefficient  is  close  to  an 
optimum  value  for  maximum  intensification.  The  result  for 
the  bulk  scheme  is  slightly  different  from  that  reported  by 
Montgomery  ef  a/.  (20 1 0) ,  where  the  decline  in  intensity  with 
increasing  drag  coefficient  occurs  at  a  much  higher  value 
than  twice  the  standard  value.  Presumably  the  difference 
here  is  because  Cd  increases  with  wind  speed  before  levelling 
off,  whereas  in  Montgomery  et  al  (2010)  its  value  was  held 
constant  for  simplicity. 
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Figure  7.  Radial  profiles  of  time-mean  azimuthally- averaged  latent  heat  flux  for  (a)  the  bulk  and  Blackadar  schemes,  and  (b)  the  MRF  and  Gayno-Seaman 
schemes  during  the  period  1 14-120  h.  Curves  1, 2, 3  denote  runs  with  Cd  =  standard/2;  Co  =  standard;  and  Cd  =  2x  standard,  respectively.  This  figure 
is  available  in  colour  online  at  wileyonlinelibrary.com/journal/qj 


The  effects  of  gradient-wind  imbalance  are  reflected  in  the 
strength  of  the  inflow,  which  in  all  schemes  is  larger  at  most 
times  for  larger  values  of  the  drag  coefficient.  However,  the 
increase  in  going  from  half  the  standard  drag  to  standard 
drag  is  typically  larger  than  for  doubling  the  standard  drag. 
In  fact,  except  with  the  Blackadar  scheme,  the  increase  with 
doubling  the  standard  drag  may  not  be  signihcant.  The 
larger  inflow  arises  because  the  increased  drag,  distributed 
through  the  depth  of  the  boundary  layer  by  subgrid-scale 
turbulent  diffusion,  leads  to  a  larger  inward  agradient  force 
in  the  layer. 

In  general,  from  an  axisymmetric  perspective,  with  a 
higher  drag  coefficient,  inflowing  rings  of  boundary-layer 
air  are  converged  farther  inwards  before  rising  out  of  the 
boundary  layer  and  ascending  into  the  eyewall  updraught. 
The  change  in  the  maximum  tangential  wind  speed  is  the  net 
effect  of  two  competing  processes:  the  increase  in  tangential 
wind  on  account  of  air  parcels  moving  to  smaller  radii 
with  partial  conservation  of  absolute  angular  momentum, 
and  the  decrease  on  account  of  the  loss  of  absolute  angular 
momentum  along  the  trajectory.  We  have  shown  that  the 
net  effect  depends  on  the  boundary-layer  scheme  used  and, 
in  particular,  on  how  diffusive  the  scheme  is  in  the  vertical 
direction.  With  the  two  schemes  that  have  a  relatively  small 
vertical  turbulent  eddy  diffusivity,  the  bulk  and  Blackadar 
schemes,  the  decrease  in  tangential  wind  in  going  from 
the  standard  value  of  drag  to  twice  its  value  is  much  less 
than  that  with  the  more  diffusive  MRF  and  Gayno-Seaman 
schemes. 

The  depth  of  appreciable  inflow  is  much  larger  with  the 
MRF  scheme  than  with  the  other  schemes.  This  result  is 
consistent  with  the  relatively  large  vertical  eddy  diffusivity 
in  this  scheme,  which  has  a  maximum  value  of  nearly 
600m^s~^  Because  the  surface  drag  is  distributed  over  a 
much  deeper  layer,  the  local  disruption  of  gradient  wind 
balance  is  smaller,  leading  to  a  substantially  weaker  inflow 
than  for  the  other  schemes. 

Maximum  values  of  vertical  eddy  diffusivity  with  the  bulk 
scheme  are  comparable  with,  but  systematically  a  little  larger 
than  those  with  the  Blackadar  scheme  and  lie  in  the  range 
60-120  m^s~^  In  contrast,  the  Gayno-Seaman  scheme  is 
relatively  diffusive  with  maximum  values  reaching  about 
250  m^  s^^  All  schemes  except  the  Gayno-Seaman  scheme 
predict  a  marked  increase  in  the  maximum  diffusivity  as  the 
drag  coefficient  is  increased.  This  increase  might  be  expected 
in  a  boundary  layer  with  mechanically-driven  turbulence  as 


in  a  hurricane.  The  lack  of  sensitivity  of  the  diffusivity  to  the 
drag  coefficient  in  the  Gayno-Seaman  scheme  is  interesting, 
but  remains  to  be  investigated. 

Despite  the  fact  that  the  exchange  coefficients  for  heat 
and  moisture  are  held  fixed,  the  finding  that  there  exists 
an  optimum  drag  coefficient  for  maximum  intensity  is 
significant  in  the  light  of  Emanuel’s  theory  for  the  potential 
intensity  of  a  tropical  cyclone  (Emanuel,  1995).  This 
theory  predicts  a  decrease  in  intensity  with  increasing  drag 
coefficient  for  a  fixed  enthalpy  exchange  coefficient.  The 
intrinsically  unbalanced  dynamics  of  the  boundary  layer 
emphasized  here  is  not  a  feature  of  Emanuel’s  theory 
(Smith  et  al,  2008;  Smith  and  Montgomery,  2008).  The 
present  study  underscores  the  importance  of  adequately 
representing  the  boundary-layer  dynamics  in  models  for 
forecasting  tropical-cyclone  intensity.  However,  it  is  not 
yet  known  which  boundary-layer  formulation  is  the  most 
appropriate  for  this  purpose,  highlighting  the  need  for  a 
concerted  research  effort  in  this  direction. 

In  summary,  we  have  isolated  what  we  believe  to  be  an 
important  aspect  of  tropical-cyclone  dynamics  and  have 
provided  a  physically  plausible  interpretation  of  the  results. 
This  interpretation  should  not  be  construed  as  a  complete 
theory.  Even  though  the  enthalpy  exchange  coefficients 
are  held  constant  in  this  study,  an  increase  in  the  surface 
wind  speed  brought  about  by  a  change  in  the  value  of 
the  drag  coefficient  results  in  a  change  in  the  enthalpy 
fluxes.  To  account  for  the  changes  in  the  aggregate  forcing 
of  the  overturning  circulation  due  to  these  flux  changes 
would  require  a  theory  incorporating  the  coupling  of  the 
(unbalanced)  boundary  layer  with  the  interior  flow.  As  far 
as  we  are  aware,  no  such  theory  exists  at  the  present  time. 
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